Optomechanical Cooling of a Macroscopic Oscillator by Homodyne Feedback 
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We propose a simple optomechanical model in which a mechanical oscillator quadrature could 
be "cooled" well below its equilibrium temperature by applying a suitable feedback to drive the 
orthogonal quadrature by means of the homodyne current of the radiation field used to probe its 
position. 

PACS numbers(s): 03.65.Bz, 42.50.Dv, 42.50.Vk 

The problem of considering a macroscopic oscillator 
in terms of Quantum Mechanics is usually avoided be- 
cause one can obtain the right results without using any 
quantum mechanical hypothesys. When, however, one 
whishes to use it as a device to detect extremely small 
displacements due to very weak forces, as in the grav- 
itational wave detectors, one has to be careful in con- 
sidering it as a mere macroscopic object. Should one 
consider a macroscopic oscillator as a quantum oscilla- 
tor, once all other possible noise sources were eliminated 
by using filters, screens, insulators etc., the ultimate cri- 
terion one has to satisfy is the one associated with the 
thermal noise [|l]j^] . For the harmonic oscillator it means 
ksT < huj m /2, where ks is the Boltzmann constant, u> m 
the mechanical angular frequency and T the temperature 
of the environment in which the oscillator lives. This 
prohibitive limit, for macroscopic massive oscillators, is 
however only valid when the measurement time r is of 
the order of the mechanical relaxation time r m . The ac- 
tual limit can be expressed as 2fcsTr/Q m < h |y]. In 
this case it is possible to consider a macroscopic mechan- 
ical oscillator as a quantum oscillator even at liquid He 
temperature ||, but very high mechanical Q m = ui m T m 
factors and also short observation times should be consid- 
ered. To have better results and, for instance, to detect 
millisecond duration bursts of gravitational waves from 
supernovae, one should measure out of resonance, as in 
VIRGO or LIGO proposals Q , or at lower temperatures, 
as in massive bar detector schemes |gj . The thermal fluc- 
tuations are, however, the fundamental limitations and, 
in order to reduce their effects, one usually should lower 
the environment's temperature. 

In this letter we present an alternative way of cool- 
ing the oscillator's observed quadrature which could be 
experimentally accessible. 

We consider an empty Fabry-Perot cavity with one 
fixed mirror with transmittivity T r and one perfectly re- 
flecting end mirror. The completely reflecting mirror can 
move, undergoing harmonic oscillations damped by the 
coupling to a thermal bath in equilibrium at temperature 
T. The cavity resonances are calculated in the absence 
of the impinging field, hence, if L is the equilibrium cav- 
ity length, the resonant frequency of the cavity will be 
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v c = lo c /2h = nc/2L, where n is an arbitrary integer 
number and c the speed of light. Furthermore, we as- 
sume that at the frequency of the impinging field vq, the 
fixed mirror does not introduce any excess noise beyond 
the input field noise. We also assume that retardation 
effects, due to the oscillating mirror in the intracavity 
field, are negligible. We shall use a field intensity such 
that the correction to the radiation pressure force, due 
to the Doppler frequency shift of the photons on the 
moving mirror, is completely negligible. This means con- 
sidering the damping coefficient of the oscillating mirror 
to be only due to the coupling with the thermal bath. 
Thus, we can write the Hamiltonian as 

H = hu c (B^B + i) + |- + \mu? m x 2 + H lnt , (1) 

where B and B^ are the boson operators of the resonant 
cavity mode; p and x are the momentum and the dis- 
placement operators, respectively, from the equilibrium 
position of the oscillating mirror with mass m and oscil- 
lation frequency v m — Lo m /2ir. The mechanical angular 
frequency Lu m will be many orders of magnitude smaller 
than oj c to ensure that the number of photons generated 
by the Casimir effect j?J is completely negligible; we actu- 
ally are in the so called adiabatic approximation, i.e., the 
cavity round trip time of the photon is much shorter than 
the mirror's period of oscillation. H int accounts for the 
interaction between the cavity mode and the oscillating 
mirror @. Since we have assumed no retardation effects, 
Hi n t simply represents the effect of the radiation pressure 
force which causes the instantaneous displacement x of 
the mirror and can be written as 

H mt = -h—B J <Bx = -hGB' l B(A + Al), (2) 
L 

where we have introduced the dimensionless ladder op- 
erators (A and A^) for the oscillating mirror, and the 
coupling constant becomes G — y/ hu> 2 / 2mu m L 2 . The 
intracavity radiation field mode B is damped through 
the output fixed mirror at a rate 7b = cT r /2L, while "f m 
is the mechanical damping rate ("f m << 7?,). 

The above interaction (||) gives rise to non lin- 
ear stochastic equations whose linearization around the 
steady state is equivalent, in a frame rotating at the im- 
pinging frequency vq, to replace Eq. (|l|) with |)| 

H 1 = hAtfb + huj m a)a + H' lnt , (3) 

where now all the operators represent small fluctua- 
tions around steady state values, i.e. B = [3 S + b 
and A = (A) s + a. These are determined by x s — 
^/h/2muj m ((A} s + (A)*) = hLo c \(3 s \ 2 /mLo 2 m L- (3 S = 
{B) s = -Jlbfiinji^i ~ with p in the classical field 
characterizing the input laser power P in = hu}Q\Pi n \ 2 . 
The parameter A represents the radiation phase shift due 
to the detuning and to the stationary displacement of the 
mirror, A = lu c — ujq — u c x s /L; further we have 
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H' lnt = h X XY v , (4) 

where X = (a + a')/2 is the mirror position quadrature, 
Y v = (be lip + b J *e~ lip )/2 is the radiation quadrature with 
tp = arg(/3 s ), and X = -4G|/3,|. Eqs. (§) and (@) repre- 
sent the starting point for further analysis of our system. 

The mirror displacement induces a phase shift on the 
radiation field, hence the latter can be used as a "me- 
ter" to measure the mirror position quadrature X. For 
simplicity, in Eq. (j^) we assume we can set A = by 
just varying the cavity detuning. This setting allows us 
to write the evolution equation for the whole density op- 
erator D as 

D = CD- 1 - [H' mt ,D] + | (2bDtf - tfbD - Dtfb) , 

(5) 

where C describes the damped dynamics of the mechan- 
ical mode a which is considered in thermal equilibrium 
at temperature T. We also assume that the number of 
thermal photons is negligibly small at optical frequency. 
Following the treatment of Refs. ]ll"||l^ ], we assume that 
the radiation mode b is heavily damped, so that due to 
the large value of 7^ >> |%|, the b mode always will be 
near to its vacuum state (we are considering fluctuations 
around the steady state (3 S ). This allows us to adiabati- 
cally eliminate the b mode and to perform a perturbative 
calculation in the small parameter x/lb fl^l - 

The measurement of the quadrature X is obtained 
by performing a homodyne measurement [(lif of a 
generic quadrature of the "meter" mode b, i.e. Y-$ = 
(be~ lS + 6^e l<5 ) /2 with 5 a phase connected to the local 
oscillator |[4| . Due to the interaction ([|) between the two 
modes, one gets information on X by directly observing 
the mode b. The continuous monitoring of the b mode 
via homodyne detection modifies the time evolution of 
the whole system. 

We are now able to apply a phase-dependent feedback 
loop to control the dynamics of the mechanical mode of 
interest a. Other feedback schemes could be devised [ jl5| , 
but only the phase-dependent feedback produces the de- 
sired effect. In order to be easily followed, we use the 
continous feedback theory recently proposed by Wiseman 
and Milburn |l6| ], who well explained the implications 
and limitations of this feedback. One has to take part 
of the stochastic output homodyne photocurrent [jl7| ob- 
tained from the continuous monitoring of the meter mode 
6, and feed it back to the mirror dynamics (for example 
as a driving term) in order to modify the evolution of the 
mode a. 

In the limiting case of a feedback delay time much 
shorter than the characteristic time of the mechanical 
mode, it is possible to obtain a Markovian equation for 
the reduced density matrix p = Tr^-D in the presence of 



feedback |16|. Thus, following Ref. fig] , we get 
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p = £p-^[X, [X, p]} + JC (ie^pX - ie-**Xp) + ^-p, 

(6) 

where T = x 2 /lb (we have defined <fi — 6 + tp, which is the 
only relevant phase of the b mode influencing the dynam- 
ics of the mechanical mode a), K. is a Liouville superoper- 
ator describing the way in which the feedback signal acts 
on the system of interest and rj represents the photodetec- 
tor efficiency. This master equation is the starting point 
of our discussion. The second term of the right hand side 
of Eq. (^|) is the usual double-commutator term associ- 
ated to the measurement of X, it results from the elim- 
ination of the radiation variables; the third term is the 
feedback term itself and the fourth term is a diffusion- like 
term, which is an unavoidable consequence of the noise 
introduced by the feedback itself. 

The mirror is considered to be in a thermal bath char- 
acterized by a damping constant 7 m , so that we have 



Cp = -iuj m [a f a, p] + 



a f , [a f - a,p] + 



l^^f[a + a\[a + a\p]] , (7) 



where [ , ]+ means the anticommutator and the limit 
ksT >> huj m is taken into account |Tq] , and due to the 
frequency we are considering, is surely valid at room tem- 
perature down to millikelvins at least. Moreover, since 
the Liouville superoperator fC can only be of Hamilto- 
nian form ]l6[ ), we choose it as fCp = g \a — a),p\ /2 
which means feeding back the measured homo- 
dyne photocurrent to the mechanical oscillator with a 
driving term in the Hamiltonian involving the mechani- 
cal quadrature orthogonal to the measured one; g is the 
feedback gain related to the practical way of realizing the 
loop. One could have chosen to feed the system with a 
generic phase-dependent quadrature, due to the homo- 
dyne current, however, it will turn out that the above 
choice gives the best and simplest result. Since the mea- 
sured quadrature of the mirror is its position the feed- 
back will act as a drive for the momentum. Using the 
above expressions in Eq. (|^) and rearranging the terms 
in an appropriate way, we finally get the following master 
equation: 

p = ~ (N + 1) (2apa f - a)ap - pa) a) (8) 
+ — N {2a) pa — aa) p — paa)) 

- ~M {2a) pa) - a) 2 p ~ pa) 2 ) 

— ^M* {2apa — a 2 p — pa 2 ) — iuj m [a^a, p] 

_ (| sm + 2^ ([ a a >p ] _ [ a v,p]) , 

where 7 = j m — g sin <f>, and 
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N = — 

7 

1 

M = — 

7 



(k B T 1\ , T , g 2 , .g 



7m U^-2j + 4 + 4^T + 2 Sm< 

* B T r 5 2 .3 
7m 7 1" 1 T~F - z o cos< 



This Eq. (Q) is very instructive because it clearly shows 
the effects of the feedback loop on the mirror mode 
a. The proposed feedback mechanism, indeed, not only 
introduces a driving term to the mirror's momentum 
quadrature, it also simulates the presence of a bath with 
nonstandard fluctuations, characterized by an effective 
damping constant 7 and by the coefficients M and N, 
which are given in terms of the feedback parameters |l2"[ . 
For the positivity of the density matrix the external pa- 
rameters should be chosen such that \M\ 2 < N(N + 1). 
This can be checked with a unitary transformation giving 
Eq. in a manifest Lindblad form |l9j for the above 
inequality. An interesting aspect of the effective bath de- 
scribed by the first four terms in the right hand side of (|^) 
is that it is characterized by phase-sensitive fluctuations, 
depending upon the experimentally adjustable phase 4>. 

Because of its linearity, the solution of Eq. (||) can 
be easily obtained, as shown in Refs. plj, by using the 
normally ordered characteristic function |2Cfl and assum- 
ing the mirror initially in a thermal state at temperature 
T, i.e. p(0) = (l- e - n »™/ kBT )Y, n \n)(n\e- nHu ™/ ksT , 
where \n) is the number state of the mode a. 

The stationary state is reached only if the parameters 
g, cj>, Lo m and 7™ satisfy the stability conditions 7 m — 
g sirup > and — 7 m <7sin0 > 0. For simplicity we 
choose <j) = — 7r/2 from now on since this choice turns 
out to be best. Under the stability conditions and in 
the long time limit (t — > 00) the variance of the position 
quadrature operator X — (a + a')/2 for the mirror is 

^ = 9 2 7m + u 2 m + j m g 



8?/r ( 7m + g){uj 2 m + i m g) 

k ^ T + r \ 7^4 ( (9) 



2hu> m 87™ / (7 m -I- g)(tJ^ + 7 m 3) 

while for the orthogonal quadrature P = (a — a')/(2i), 
i.e. the mirror's momentum, we get 



(P 2 



g 2 u 2 m 



8r)T ( 7m +g){uJ 2 a + J m g) 
( k B T r \ g 2 +<+ 7m g 

— : 7 m — — 



+ h^ + 7— 7m , " ' " m 2 r y (10) 

In the case of no coupling with the cavity mode, the 
above variances for the macroscopic oscillator only con- 
sist in the thermal noise as one should expect. When- 
ever an indirect detection of the mirror position is made 
the backaction noise is added. The latter, however, is 
usually negligibly small compared with the previous one. 
Instead, by using the feedback with sufficient high gain 
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(g >> oj m Q m ), we can set T e ff w Tuj^/g 2 as an effec- 
tive temperature, and the mirror's position quadrature 
variance becomes 

&y m g* 8^ * ' 

Although the feedback introduces excess noise, it also 
gives a scale factor for the thermal noise term by means 
of T e ff, so that with an appropriate choice of the param- 
eters, the latter can be strongly reduced. 

It is also to remark that the proposed phase-dependent 
feedback does not produce a proper squeezing; moreover, 
it can extract the thermal noise from the system, because 
the variance reduction occurs, for not extremely high val- 
ues of g, in both quadratures, as can be evicted from Eqs. 
(|| |l(]). Hence, it acts as a refrigerator. 

To better show the potentiality of this feedback mecha- 
nism let us consider the spectrum of the position quadra- 
ture. To this end the Fourier transforms of the stochas- 
tic equations connected with the master equation (p|) are 



easily written down 21 



iu)X(u) = w m P{w) - gX(uj) - ^jX in (uj) 

iu)P(u) = ~uo m X(uo) - 7 m P(w) - ^fP in {oj) (12) 

and the input noise operators have the following correla- 
tions 

(X m (cj)X in (-uj r )) = i(2JV + 1 + 2Re{M})<5(w - u') 
(P m (uj)P m (-cj r )) = i(2JV + 1 - 2Re{M})S(cj - J) 
{X in {u)P in (-J)) = ~(i + 2Im{M})<5(w - w') . (13) 

Defining S g {u>) = {X(lu)X(—uj)) s , we get from Eqs. (|l 

( 7r 2 „+c 2 +^ 2 „)(2A + l) 



5 sH = 77^ 1 



4 

( 7 l+u 2 -ul)2Re{M} 



(14) 



where the subscript s indicates the symmetrized correla- 
tion and S(w) = [(iu> + g)(iu> + J m ) + w 2 n ] . 

As a practical example we take the physical parameters 
of the model presented in Ref . JT^] ■ Taking their values 
for granted, i.e. m = 10 Kg, v rn = 10 Hz, j. m = 1 s _1 , 
L = 4 m, v = 5.82 x 10 14 Hz, T r = 0.02, P m = 10 W, we 
get r w 200 s _1 , then x ~ 10 4 s" 1 . This choice satisfies 
the relation "ft >> X an d au others inequalities, and we 
further take rj w 1 and T = 300 A". 

Then, in Fig. 1 we show the (scaled) spectrum of Eq. 
( |j"4| ) for various values of g. The curve for g — practi- 
cally coincides with the spectrum of the mirror not cou- 
pled to the cavity mode because of the smallness of the 
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backaction noise. It results evident that for high values 
of the feedback gain the spectrum is practically vanish- 
ing while the peak at the mechanical resonance frequency 
gradually disappears and one peak at zero frequency ap- 
pears, with very small amplitude. With the proposed 
feedback, a transition from a dissipative to a diffusive be- 
haviour of the oscillator is obtained by just varying the 
feedback gain. These results, although in a different con- 
text, are similar, but not equivalent, to those obtained in 
Rcf . p| , where the feedback was used for the regulation 
of a microcantilever response and a direct photodetec- 
tion was used, instead of our phase-dependent scheme. 
Furthermore, the temperature of the bath was assumed 
negligibly small. 

Summarizing, we have proposed a feedback scheme 
based on an indirect measurement to reduce position 
quadrature uncertainty of a macroscopic oscillator. The 
described mechanism can be very useful in reducing the 
effect of thermal noise in quadratures of macroscopic mir- 
rors, as those devoted to the gravitational wave detectors, 
even at room temperatures. The feedback loop may con- 
sist in a transducer Q which transforms the random op- 
tical signal in a stochastic electric signal which in turn 
acts as a mechanical driving on the mirror's momentum. 
This could be readly realized, but it is not the only way, 
by using the feedback current to vary the potential of 
a capacitor formed by the oscillating mirror and a fixed 
plate. On the other hand, depending on the specific ex- 
perimental realization of the feedback loop there could 
be some limitations on the values of g. 

We think that the practical implementation of the dis- 
cussed model, even though in a situation far from the 
oversimplified theoretical one, should be an interesting 
challenge for an experimentalist, and it will turn out ex- 
tremely useful in reducing the thermal fluctuations with- 
out lowering the bath temperature. 
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FIGURE CAPTIONS 

Fig. 1 The quantity S g (u) / '(2n(X 2 ) g=0 ) is plotted (in 
a semilog scale) versus w for the following values of g in 
s" 1 : a) 0; b) 1; c) 10; d) 10 2 ; c) 10 3 . We have used the 
values of other parameters listed in the text. 
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